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O ! Abstract 



Q^ . We determine the sensitivity of the KamLAND and Borexino ex- 

^^ [ periments to the neutrino regeneration effect in the Earth as a function 

1-^ ' of Am? and 9, using reahstic numbers for the signal and background 

']P' rates. We compare the results obtained with the x^ method with those 

Qh. obtained from the conventional day- night asymmetry analysis. We 

(— I I also investigate how well one should be able to measure the neutrino 

oscillation parameters if a large day-night asymmetry is observed, 
taking the LOW solution as an example. We present an enlarged 
rS I parameter space, which contains mixing angles greater than 7r/4 where 

C^ ' the heavy mass eigenstate is predominantly fe, and determine the 

electron neutrino survival probability for this traditionally neglected 
scenario. We emphasize that this portion of the parameter space yields 
different physics results when dealing with the MSW solutions to the 
solar neutrino puzzle and should not be neglected. 
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1 Introduction 

A number of experiments [l|, 0, |^, ||, |^ have accumulated over the years 
a large amount of solar neutrino data. The data indicate that the num- 
ber of solar neutrino induced events is significantly smaller than expected 
and, furthermore, that the electron neutrino survival probability is energy 
dependent. This "solar neutrino puzzle" is best solved by assuming that 
the electron-type neutrino [ug) oscillates into another active neutrino species 
(some linear combination of the muon-type neutrino (z/^) and the tau-type 
neutrino (z^r))? or a sterile (weak isosinglet) neutrino. In light of the very ro- 
bust Super-Kamiokande evidence for z/^ atmospheric neutrino oscillations Q , 
the oscillation of solar neutrinos seems a very likely and natural hypothesis. 

The current experimental situation is such that there are four discon- 
nected regions in the two-neutrino oscillation parameter space that fit the 
data. One of them, the "just-so" solution, relies on vacuum neutrino oscil- 
lations with a very long wavelength (comparable to the Earth-Sun distance) 
0, while the other three 0, |^ rely on the MSW effect [^] to produce the 
required energy dependence of the electron neutrino survival probability. 
Discriminating among all these solutions is the goal of the current and the 
next generations of neutrino experiments. 

Even though one can classify the solar neutrino puzzle as strong evidence 
for neutrino oscillations, it is as yet not considered definitive. The main foci 
of criticism traditionally have been that the Standard Solar Model (SSM) 



101 might not be accurate enough to precisely predict the fluxes of different 
energy components of solar neutrinos, and that the evidence for solar neutrino 
oscillations relies on a combination of hard, different experiments. Even 
though it seems very unlikely that reasonable modifications to the SSM alone 



can explain the current solar neutrino data (see, for example, ||TT|), one still 
cannot completely discount the possibility that a combination of unknown 
systematic errors in some of the experiments and certain modifications to the 
SSM could conspire to yield the observed data. To conclusively demonstrate 
that there is indeed new physics in solar neutrinos, the experiments now are 
aiming at detecting "smoking gun" signatures of neutrino oscillations, such as 
an anomalous seasonal variation in the observed neutrino flux or a day-night 
variation due to the regeneration of electron neutrinos in the Earth. In this 
paper we study the sensitivity reach of two upcoming neutrino experiments, 
Borexino and KamLAND, to the Earth regeneration effect. 

Out of all solar neutrino components, both experiments will be most 



sensitive to ^Be neutrinos. These are neutrinos produced in the electron 
capture by '^Be nuclei in the Sun's core ("^Be+e" — >'^Li +z/e). One very 
important fact is that these neutrinos are almost monochromatic, with E,^ = 
0.862 MeV (90% of the time) or E^ = 0.383 MeV (10% of the time), 
depending on the final state of the ^Li nucleus. Since the Ei, = 0.383 MeV 
neutrinos cannot be cleanly seen in future detectors, we will only consider the 
El, = 0.862 MeV neutrinos, which will be referred to as the ^Be neutrinos. 

The study of the ^Be neutrino flux is particularly important, for a variety 
of reasons. First, in the SSM independent analysis of the solar neutrino data 



rT| , where one arbitrarily rescales the flux of neutrinos from different sources, 
the flux of ^Be neutrinos comes out extremely suppressed (in fact the best fit 
value for the ^Be flux is negative!), and the measurement of a reasonable flux 
would dramatically constrain such attempts. Second, since the prediction 
of one particular MSW solution (the small angle solution) for the survival 
probability of ^Be neutrinos is very different from the other two solutions, 
one can separate if from the other two by measuring the ^Be solar neutrino 
flux. Third, as was recently shown |]T2[, one can either establish or exclude 
the "just-so" solution by analyzing the seasonal variation of the ^Be solar 
neutrino flux at Borexino or KamLAND. Finally, it might also be possible 
to separate the Ue from the u^^r component in the ^Be flux, by studying the 



kinetic energy spectrum of recoil electrons |]T3 in future experiments. 



It has been known for over a decade that the propagation of solar neutri- 
nos through the Earth can result in a measurable variation in the observed 
neutrino event rates |]14|. Reference |]15|, in particular, contains a detailed 



analysis of the expected day-night asymmetry for the Super-Kamiokande, 
Borexino, and SNO experiments. In this paper we extend the previous anal- 
yses in several important aspects. First, we present an enlarged parameter 
space, where the vacuum mixing angle is allowed to vary over its entire 
physical range from to 7r/2Q We flnd that not only does the day-night 
asymmetry stay nonvanishing at maximal mixing {6 = vr/4), in agreement 



with ^7j, but that it also smoothly extends into the other part of the 
parameter space (7r/4 < 9 < 7r/2). Second, we display the sensitivity regions 
of KamLAND and Borexino in this enlarged parameter space, using realistic 
numbers for the signal and background rates. In our analysis we use the x^ 
method, and study the effect of various binning schemes. Finally, we explore 



*This enlarged parameter space has already been mentioned in the context of three- 
flavor oscillations ilql. 



the possibility of using the neutrino regeneration data at the two experiments 
in question to measure the oscillation parameters. 

This paper is organized as follows. In Sec. ^ we review the day-night 
effect and present the day-night asymmetry expected for ^Be neutrinos as a 
function of the two neutrino oscillation parameter space. We also introduce 
an enlarged parameter space, < ^ < 7r/2. In Sec. |^, we study the sensitivity 
of the KamLAND and Borexino experiments to the day-night asymmetry 
and to the zenith angle dependence of the ^Be flux. In Sec. ^ we study the 
possibility of measuring the oscillation parameters if a significant day-night 
effect is observed at either Borexino or KamLAND. We contrast the analysis 
of the day-night asymmetry with the zenith angle distribution. In Sec. |^ we 
present a summary of our results and conclusions. 

2 Electron Neutrino Regeneration in the Earth 

As was realized over a decade ago 0, neutrino-matter interactions can dra- 
matically affect the pattern of neutrino oscillations. The reason for this is 
that neutrino- matter interactions are flavor dependent, given that the matter 
distributions of interest (the Earth, the Sun) contain only first generation 
particles. One well-known consequence of this is that, in the case of neutrinos 
produced in the Sun's core, it is possible to obtain an almost complete 
^e ~^ ^'other transformation even when the vacuum mixing angle is very small 



It has also been pointed out by several authors ||I^ |T5|, |17| that matter 
effects might also be relevant for neutrinos traversing the Earth. One ex- 
perimental consequence of neutrino-Earth interactions is that the number of 
events detected during the day (when there are no neutrino-Earth interac- 
tions) can be statistically different from the number of events detected during 
the night. The Super-Kamiokande experiment has already presented exper- 
imental data which seem to slightly prefer a nonzero day-night asymmetry, 
even though the result is not yet statistically significant |]TB[ (the most recent 
result is Ajd^ = 0.065 ± 0.031 ± 0.013). 

In this section we review the electron neutrino regeneration effect in 
the Earth and how it modifies the solar neutrino data. We also present 
the expected day-night asymmetry for ^Be neutrinos at the KamLAND and 
Borexino sites. 



2.1 The Day-Night Effect 

If neutrinos have mass, it is very likely that, similar to what happens in 
the quark sector, neutrino mass eigenstates are different from neutrino weak 
eigenstates. Assuming that only two neutrino states mix, the relation be- 
tween mass eigenstates and flavor eigenstates is simply given by 

\ui) = cos6'|z/e) — sin6'|z/^), 

lz/2) = sin6'|i^e) + cos^|z/^), (2.1) 

where 6 is the vacuum mixing angle, {ui) and |z/2) are the mass eigenstates 
with masses mi and m2, respectively, and z/g *-> u^ mixing is considered. The 
mass-squared difference is defined as Arn^ = m^ — m\. We are interested in 
the range of parameters that encompasses all physically different situations. 
First, observe that Eq. ( p.l| ) is invariant under 9 ^ 9 + Ti, u^ ^ —I'e, ^fi ^ 
—Vixi i-G. 9 E [— 7r/2,7r/2] and 9 G [7r/2,37r/2] are physically equivalent. 
Next, note that it is also invariant under 9 -^ —9, u^ -^ —v^i-, ^2 —>■ —^2, 
hence it is sufficient to only consider 9 G [0, 7r/2]. Finally, it can also be made 
invariant under 9 -^ 7r/2 — 9,1/^-^ —v^ by relabeling the mass eigenstates 
Vi •«->■ z/2, i.e. Am"^ —>■ —Am?. Thus, all physically different situations are 
obtained if < sin^ 9 <1 and Am? is positive, or < sin^ 9 < 1/2 and Am^ 
can have either sign. In what follows, we will use the first parametrization 
(Am^ > 0), unless otherwise noted. 

^Be neutrinos reach the Earth as an incoherent mixture of |i/i) and 11/2) 
| T9| (see also pO|, |17|); with probabilities Pi and P2 = 1 — -Pi as long as Am^ > 



10^^ eV^. Pi is given in Eq. ( |A.7| ) in terms of the jumping probability Pc and 
its value depends on the details of the neutrino production and propagation 
inside the Sun, as presented in Appendix 0. The probability Pee of detecting 
a Ue on the Earth is given by 

Pee = PlPle + (1 " Pl)P2e , (2.2) 

where Pie is the probability that vi (z/2) is detected as a Ve for i = 1 (2). 
Because Pie + P2e = 1 (always, independent of matter effects, because of the 
unitarity of the Hamiltonian) , one can rewrite Eq. ( p.2|) 



Pee = Pi + (1 - 2Pi)P2e . (2.3) 

In the case of neutrinos detected during the day, P2e = sin^ 9 (the vacuum 
result), while for neutrinos that traverse the Earth P2e = Pg^ must be 



calculated numerically, and depends on the density profile of the Earth 
and the latitude of the location where the neutrinos are to be detected. 
One should also remember that muon or tau neutrinos still interact in the 
detector through neutral currents, although the even rate is down by a factor 
of i? ~ 0.2 compared to electron neutrinos. The day-night asymmetry 
{Adn = (events detected during the night minus events detected during 
the day) /(total)) is, therefore, 

^^ (2Pi + (1 - 2Pi)(P2e + sin2 0))(l ~R) + 2R- ^ ' ' 

It is important to note that A^n does not have to vanish, as used to be 
the general lore in the past, when sin^ ^ = 1/2 (maximum mixing), as was 
clearly shown in |]T^. A^jsr does vanish, of course, when Pi = 1/2 (a fifty-fifty 



mixture of mass eigenstates reaches the Earth). 

It is interesting to note that, in the past, Adn was always computed 
assuming that sin^ 6^ < 1/2. However, it is perfectly acceptable to have 
sin^ 9 > 1/2, when the heavy mass eigenstate (z/2) is predominantly Ue- 
While in the case of vacuum oscillations physical results depend only on 
sin^ 29, in the case of neutrino-matter interactions sin^ 9 > 1/2 leads to 
physically different results. Using sin^ 29 as a parameter in the latter case 
can be misleading, as < sin^ 2^ < 1 does not cover all physically distinct 
possibilities. Similar to what was pointed out in [^ for the transition 
between sin^ 9 < 1/2 to sin^ 9 = 1/2, we will show that for the entire range 
of < sin^ 6' < 1 the behavior of A^^ is smooth. In Appendix ^ we explain 
in detail how to extend the expression for Pi to the case sin^ 9 > 1/2. 

2.2 The Day-Night Asymmetry at 36° and 42° North 

We numerically compute the value of Pg^ and A^jy for ^Be neutrinos at 
KamLAND (latitude = 36.4° north) and Borexino (latitude = 42.4° north). 
We assume a radially symmetric exponential profile for the electron number 
density inside the Sun, and use the analytic expression for the survival prob- 



ability of neutrinos produced in the Sun's core derived in [^, as presented 
in Appendix ^. We appropriately integrate over the ^Be neutrino produc- 
tion region inside the Sun, using the results of the SSM [|1^, conveniently 
tabulated in |2^. 

We use a radially symmetric profile for the Earth's electron number den- 
sity, given in 12^ , and the zenith angle exposure function for the appropriate 



latitude, which was obtained from |^^. For a plot of the electron number 
density profile in the Earth see Fig. 2 in [|1^ and for the zenith angle exposure 



function see the upper left-hand corner of Fig. 5 in |T^ . The model predicts 
that the electron number density in the Earth's mantle varies in the range 
2.1 to 2.7 moles/cm^, while in the outer core the electron number density 
is significantly greater (4.6 to 5.6 moles/cm^). Because of the latitude of 
Borexino and KamLAND, the solar neutrinos detected at these experiments 
will not travel through the inner core. 

Fig. m depicts the constant day-night asymmetry contours for '^Be neu- 
trinos]] at KamLAND and Borexino. It is important to note that, unlike 
conventionally done in the literature, the x-axis here is sin^ 0, not sin^ 26. To 
facilitate comparison with earlier results, we also depict the same information 
in the (Am^,sin^2^) plane in Fig. |^, where once again we vary the mixing 
angle in its entire physical range < ^ < 7r/2. 

As Fig. 1^ demonstrates, the asymmetry contours smoothly extend into 
the sin^ 6' > 0.5 half of the parameter space. One can see that in that region 
the day-night asymmetry is non-zero and may, in fact, be quite large. This 
kind of behavior had already been seen in |l^, for day-night asymmetry 
contours at SuperKamiokande (see Fig. 11 in [1^). This is to be contrasted 
with conventional analyses, which choose axes as in Fig. 0, but only show the 
< ^ < 7r/4 half of the parameter space. As a result, contours there seem 
to abruptly terminate at maximal mixing. 

It is also easy to see from our plots that, with the choice of variables 
as in Fig. |l], there is nothing special about maximal mixing. This point is 
somewhat obscured in the (Am^, sin^ 26) plane, where it seems that the slope 
of the contours abruptly changes around sin^ 2^ = 1. The reason for this is 



that the Jacobian of the transformation from sin^ 6 to sin^ 26', 

d[sui26) ^ 2^^^ (2 5) 

d{sni6) cos 6 

vanishes at maximal mixing 6 = 7r/4. It can be argued, therefore, that sin^ 6 
represents a more natural parametrization. From here on we will always use 
sin^ 6* as a parameter .0 

^We only assume Ve oscillations into other active neutrino species. 

■flf one wishes to keep the symmetry between 9 < 7r/4 and 9 > 7r/4 for vacuum 
oscillations while avoiding the singular Jacobian, the best choice for the horizontal axis 
would be tan^ in log scale, as was done in |l6| in the context of three-flavor oscillations. 
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Figure 1: Constant day-night asymmetry contours (10%, 5%, 1%, 0.5%) in 
the (sin^ 9, Am^)-plane for ^Be neutrinos at the KamLAND and Borexino 
sites. The vertical dashed hne indicates sin^ 9 = 1/2, where the neutrino 
vacuum mixing is maximal. 
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Figure 2: Constant day-night asymmetry contours (10%, 5%, 1%, 0.5%) in 
the (sin^ 29, Ar7i^)-plane for ^Be neutrinos at the KamLAND and Borexino 
sites. The right side of the plot, with decreasing scale, can also be thought 
of as Am^ < 0, 9 < 45°. 



The day-night asyminetry for 9 = it /A is in general non-zero and, indeed, 
can be larger than 10%. Our analysis, thus, is in complete agreement with the 



findings of ^7j and extends them to the other half of the parameter space. 



Note that constant day-night asymmetry contours do close as sin 9 —> 1. 
This is expected, because in that limit, just like for sin^ 9-^0, there is no 
neutrino mixing, and so Pee goes trivially to 1 and A/jjv vanishes. 

Almost all other features of the contours in Figs. |I] and ^ can also be 
understood analytically. Several physical effects are involved in shaping up 
the contours. In the low Am^ region the oscillation length in the Earth is 
comparable to the size of the Earth, independent of the value of Arn^. This 
can be understood very easily in the approximation that the Earth's electron 
density is uniform. In that case the neutrino oscillation length is given by 



n 



^r {V2GFN,y - 2 



2E, 



\2E, 



V2GfN^cos29 



or numerically 



-1/2 



(2.6) 




N. 



1 mole/cm^ 

n-l/2 



COS 2^ 



(2.7) 



and, for AmV(2E^) < V^GpN^, Lose ^ 8.2 ■ 10^ km x (1 m.o\e/cm^ /N^). 

For very small Am^ the asymmetry vanishes for two reasons. First, 
the MSW transition inside the Sun becomes non-adiabatic. For Am? -C 
10-5 eV2, 9m ^ vr/2 (Eq. (Q) in the Sun's core and Pi ~ P^ (Eq. (^)). 
As the value of the jumping probability Pc changes from to cos^ 9 it passes 
through 1/2 (for 9 < 7r/2) where A^jq vanishes, according to Eq. ( p.4| ). As 
can be deduced from Eq. (|A.16|) , the contours of constant jumping probability 
Pc are approximately described by Aw? sin^ 9 = constant, provided sin^ 9 <C 
1 and Am^ ^ 10"^ eV^. Second, the mixing angle in the Earth becomes 
close to 7r/2 and no regeneration takes place in that limit (see also Eq. ( |2.8|) 
below, where 9m — > 7r/2 gives P2e -^ sin^^). Below the line Pi = 1/2 the 
asymmetry is negative and very small. 

In the region Arra^ ^ 3 x 10^^ eV^ neutrinos undergo many oscillations 
inside the Earth, as can be seen from Eq. ( p.7|) . The relevant quantity in 



this case is the average survival probabihty, obtained after integrating over 
the zenith angle. One can understand the shape of the asymmetry contours 
in this region by, once again, approximating the electron number density in 
the Earth by a constant value. In this model, it is easy to show that, if a 
state |z/j) enters from vacuum into the Earth, the average survival probability 
inside the Earth is 

P2T = sin' Om + sin' {9 - 9m) cos 29m, W^) = \ ^2) , .. .^ 

Pf; = sin' 9m + cos'(^ - 9m) cos29m, Wi) = Wi). ^ ^ 

Here 9 is the mixing angle in vacuum and 9m is the mixing angle inside the 
Earth (see Eq. (|A.5|) ). Obviously, P^^ + P2e = 1. Using these expressions, 
one can compute the day-night asymmetry for this simplified model: 

^ ^ PN + il-PN)R-PD-il-PD)R , . 

Pn + {I - Pn)R + Pd + {I - Pd)R' ^ ■ ' 



where 



•d = sin'^0((l-P,)sin'0 + PcCos'0) + 
cos' ^o((l - PJ cos' 9 + P^ sin' 9), 
: sin'^o((l-P,)P- + P,Pf:) + 

cos' ^o((l -P,)P- + P,P-)- 



N 



9q denotes the mixing angle at the production region in the core of the 
Sun, Pc the jumping probability (Eq. ( |A.16D ), and P is a contribution of 
v^T- interacting through the neutral current interactions in the detector. We 
found that for A^e ~ 3 — 4 moles/cm^ the contours of constant A are in 
good agreement with the day-night asymmetry contours in Fig. |^ for Am' ^ 
3 X 10-*^ eV'. 

Using this simple model we can explain the behavior of the asymmetry 
contours in the large Am? region. For example, according to Fig. |l], as 
sin' 6' decreases for fixed Am', the value of the asymmetry goes down. This 
happens because, while the difference in the numerator of Eq. (|2.9| ) goes to 
zero, the denominator approaches a constant value due to the non-vanishing 
neutral current contribution. Notice that in a real experiment, in addition 
to the neutral current contribution, there will be a term proportional to 
the rate of background events, further decreasing the sensitivity. Thus, 
using asymmetry contours in this region to read off the sensitivity can be 

10 



misleading. This would be even more obvious in the case of oscillations to a 
sterile neutrino. We will return to this issue in the next section. 

Even more subtle features can be understood within this model. For 
instance, we found that the slight change of the slope seen for the 0.5% 
contour around sin^6' ~ 0.04 is due to the significant deviation of the value 
of Oq from 7r/2 in that region. 

Finally, in the region Am^ ~ 10~^ eV^ the regeneration efficiency ex- 
hibits a very strong zenith angle dependence. Because the magnitudes of 
/S.m? / {2Eu) and \/2GFNe in the core are almost equal, the mixing in the 
core is almost maximal {6m ~ 7r/4, see Eq. ( |A.5| )), while in the mantle it 
is small {6m ~ 7r/2). As a result, for neutrinos traveling through the outer 
core the conversion into u^ is much more efficient than for ones going only 
through the mantle. The oscillations do not average out completely in this 
case, resulting in the presence of several wiggles. We have explicitly checked 
that these wiggles are not washed out by the effect of the finite width of the 



^Be line l24 



Our results for 6 < tt/A agree qualitatively with the results presented in 
[ll5| for the Borexino site. The agreement is not complete, however. For 



instance, the contours in |T^ do not exhibit any wiggles in the range Am^ 
10-6 eV^. 



3 The Neutrino Regeneration Effect at Kam- 
LAND and Borexino 

In this section we study the sensitivity of the KamLAND and Borexino 
experiments to the day- night effect. 

Borexino [^ is a dedicated ^Be solar neutrino experiment. It is a large 



sphere containing ultrapure organic liquid scintillator (300 t) and can detect 
the light emitted by recoil electrons produced by elastic u-e scattering. By 
looking in the appropriate recoil electron kinetic energy window, it is possible 
to extract a very clean sample of events induced by ^Be neutrinos, if the 
number of background events is sufficiently small. Borexino expects, in the 
absence of neutrino oscillations, 53 neutrino induced events/day according 
to the SSM, and 19 events/day induced by background (mainly radioactive 



impurities in the detector, see O, E5| for details) 



The KamLAND experiment, located in the site of the original Kamiokande 



11 



experiment, was initially designed as a reactor neutrino experiment. Re- 
cently, however, the fact that KamLAND might be used as a solar neutrino 
experiment has become a plausible and exciting possibility [^ . 

KamLAND is also a very large sphere containing ultrapure liquid scintil- 
lator (1 kt), and functions exactly like Borexino. The outstanding issue to 
determine if KamLAND will study solar neutrinos is if the background rates 
can be appropriately reduced. KamLAND expects, in the absence of neutrino 
oscillations, 466 neutrino induced events/kt/day according to the SSM, and 
217 events/kt/day induced by background (mainly radioactive impurities in 
the detector, see [0, ^ for details). We will consider a fiducial volume 
of 600 t, so that 280 (unoscillated) signal events/day and 130 background 
events/day are expected. We assume that the number of background events 
is constant in time. 

We generate a histogram of the number of events expected in each of the 
A^ day and N night bins for different values of (Am^, sin^ 9). The number of 
events per year in the i-th bin is 

/events\ /days\ ,, / ,■ / ,• n nx /events\ „ 

"• i^sr) = ''' {^) (''-+^-(n-.+(i-i'i)fl)) [-^) /. (3.1) 

where Srate = 280 (53) events/day and 6rate = 130 (19) events/day for 

KamLAND (Borexino), P^^ is the electron neutrino survival probability in 

the i-th bin, R is the ratio of the i/g-e to z/^,r-e elastic cross sections^ (see 

g, at KamLAND (Borexino) R = 0.214 (0.213)) and /i=(size i-th bin 



divided by the sum of the sizes of all the bins), such that X)?^ /« = 1- As an 
example, if there are 24 (12 day, 12 night) hour-bins, f] = 1/24 for all i. In 
reality, we are interested in zenith angle bins, and in order to the determine 



fi, the exposure function presented in |T5[ is used. Note that we assume only 
statistical uncertainties. 



X 



^ is defined as 



N /"flight _ ^day^^ 



/ nig 

^' = E , A^,2 '; ^^ ^2 + N. (3.2) 



^=l ( . /^night ^ /day 



The factor A^ is included in the definition of x^ i^ order to take statistical 
fluctuations of the data into account. A detailed explanation of the philoso- 



phy behind this procedure can be found in [12 



"In the case of Ve <-> ^'sterile oscillations, R = 0. 

12 



It is important to comment at this point that, in hght of the definition of 
X^ (Eq. ( p.2|) ), the sensitivity of the experiments to the Earth matter effect 



does not require any input from the SSM, including the ^Be solar neutrino 
fiux, or from a direct measurement of the background rate. This is because 
we are comparing the night data with the day data, and no other inputs are 
required. Our quantitative results, however, depend on the expected number 
of signal and background induced events, since these quantities are used as 
input for the "data" sample. 

We will define the sensitivity of a given experiment to the Earth matter 
effect by the value of x^, computed according to Eq. (|3.2|). The sensitivity 
defined in this way depends clearly on A^, the number of day and night bins, 
and on /j (see Eq. ( p. ID ), or on the "size of the bin" . With the real experimen- 
tal data, one will certainly consider many different types of analyses in order 
to maximize the sensitivity of the data to the neutrino regeneration in the 
Earth (options include computing moments of the zenith angle distribution, 
Fourier decomposing the data, maximum likelihood analysis, and others), 
but, since we analyze thousands of "data samples" (one for each value of 
(Am^,sin^ ^)), this simple x^ approach will suffice. 

We consider two options for the size of zenith angle bins. In one of them, 
each bin has the same size, that is, the bins are equally spaced {e.g. 0° — 
30°, 30° — 60°, 60° — 90°, etc). The other option is to choose the bin size such 
that the distribution of the day data is uniform. It is worthwhile to comment 
that the latter scheme may be considered the most natural one for KamLAND 
and Borexino, which are real time experiments with no directional capability. 
In these experiments, it is straightforward to organize the data into time bins, 
which then have to be translated into zenith angle bins by associating the 
time of the event with the position of the Sun in the sky. 

Another issue to consider is the value of A^ which optimizes the sensitivity. 
It is clear that for A^ = 1 (the day-night asymmetry case) the statistical 
significance is enhanced for overall changes in the number of events, but for 
larger A^, one should be more sensitive to distortions in the zenith angle 
distribution. Different binning schemes of the "data" for Am? = 1.12 x 
10~^ eV^, sin^ 9 = 0.398 and three years of KamLAND running are depicted 
in Fig. ^, for A^ = 1, A^ = 10 equally spaced bins, and A^ = 10 "uniform" 
bins.[| 

TThe residual non- uniformity seen in the figure is due to the fact that we used a discrete 
table of values for the exposure function. 



13 



xlO' 




1000 



20 40 60 80 100 120 140 160 180 



40000 
30000 
20000 
10000 




-i-i-p-i-i-i-p-i-i-i- 

Error = 1 o times 10 



-I— T— T— T— I— r- 



hH 



* 



.. I ... I ... I ... I . I . I ... I ... I ... I 



20 40 60 80 100 120 140 160 180 



25000 
20000 
15000 
10000 



-I— I— I— I— I— i- 



Error = 1 a times 5 



** 



* 



I — i- 



-i-i-i-i-*%***t*+*|* 



•Ift-H 



-i-^*^-i- 











20 40 60 80 100 120 140 160 180 



Figure 3: Different binning schemes, for Am^ = 1.12 x 10^^ eV^, sin^ ^ = 
0.398: (a) A^ = 1 bin (the day-night asymmetry), (b) A^ = 10 equally spaced 
zenith angle bins, and (c) A^ = 10 "uniform" bins, where the day-time 
data is (roughly) uniformly distributed. The error bars contain statistical 
uncertainties only. We assume three years of KamLAND running. 
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Fig. ^ shows a comparison of the sensitivity reach of KamLAND after 
three years of running for two different binning schemes, A^ = 1 vs. A^ = 10 
"uniform" bins. The contours are drawn at 95% CL. One can easily see 
that for most of the parameter space, the best sensitivity is reached with 
the A^ = 1 case, while for a small region in the parameter space, when 
sin^ 6' ^ 0.1 and lS.rv? ~ 10~® eV^, the A^ = 10 scheme is more successful. 



This result is consistent with the analysis of Section |2.2| . As explained there, 
for IS.W? ~ 10~^ eV^ the data shows a large enhancement in the low zenith 
angle bin, while little effect in other bins. At Borexino this effect will be 
somewhat less pronounced because it is farther from the Equator. 

One can see that the contours in Fig. are similar in shape to the 



day-night asymmetry contours of Section ^J\ but quantitatively different. 



One important difference is that for !S.vn? ^ 10~^ eV^ the x^ contours 
do not extend as far in the low sin^ Q region as the asymmetry contours. 
While for low Am^ the 95% CL. contour corresponds to the day-night 
asymmetry of roughly 0.5%, for Am^ ^ 10^^ eV^ the corresponding value 
of the day-night asymmetry is at least two times greater. This phenomenon 



was already mentioned in Section 2.2. The difference occurs because the x 



2 



analysis includes, in addition to the neutral current interactions, the constant 
background rate, thus eliminating the major shortcoming of the day-night 
asymmetry analysis. 

In order to present the final sensitivity reach of KamLAND and Borexino, 
we combine the confidence level contour obtained in the different types of 
analyses, with different number of bins. Fig. ^ depicts the "optimal" 95%, 
3cr, and 5(T confidence level (CL.) contours for the sensitivity of three years 
of KamLAND and Borexino data to the day-night effect. The confidence 
levels are optimized by considering the union of same CL. contours for all 
values of A^ and both binning schemes. The day-night asymmetry provides 
the best sensitivity reach for most of the parameter space, while the A^ ?» 10 
uniform bins scheme at KamLAND increases the sensitivity for particular 
regions of the parameter space, as was discussed earlier. 

Fig. ^ clearly demonstrates that in the case of the LOW MSW solution 
to the solar neutrino puzzle, both KamLAND and Borexino should be able 
to see a larger than 5(T effect, while in the case of the SMA no significant 
effect should be detected.^] Both experiments are sensitive to a large portion 
of the parameter space which extends into Q > 7r/4 region, where the heavy 

■^KamLAND may also be sensitive to a very small portion of the LMA solution. 
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Figure 4: Comparison of the sensitivity reach of three years of KamLAND 
running with 1 bin and 10 uniform bins. 

neutrino eigenstate is predominantly z/g. 

On the other hand, should no regeneration effect be observed, a large 
portion of the parameter space, including the entire LOW region might be 
excluded. The exclusion will require knowledge of the ^Be neutrino flux, 
which can be measured, for example, by studying the seasonal variation of 
the observed event rate |]T2[. If the flux measured in this way turns out large 
and no day-night asymmetry is observed, one will be able to exclude the 
LOW solution without relying on the solar model. If, however, the measured 
flux is very small, the exclusion will be solar model dependent. 

Since the sensitivity of Borexino (KamLAND) to the day-night asymme- 
try goes down to the 1.5% (0.5%) level, it is important to consider systematic 
effects in this measurement. It is, however, difficult to anticipate systematic 
uncertainties in the absence of data. We instead looked at the measurement 
of day-night asymmetry at SuperKamiokande [|l^. The dominant systematic 
uncertainty there is the possible asymmetry in the detector, giving ±0.6%.| 
Because the recoil electrons from ®B neutrinos are forward peaked, the day 
(night) time data are detected primarily by the lower (upper) half the detec- 



§Note that the talk m ||T^ Hsts the systematic uncertainties in D/N ratio, which are 
twice as large as uncertainties in the asymmetry {D — N)/{D + N) « {{D/N) — l)/2. 
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Figure 5: 95% (darkest), 3a (dark), and 5cr (light) sensitivity confidence level 
(CL.) contours for three years of KamLAND running. The LOW solution, 
which extends from Am^ ~ 3 x 10~® eV^ to Am^ ~ 3 x 10"'^ eV^ and has 
sin^ 9 ~ 0.3 — 0.5 M, is completely covered at more than 5cr CL. 
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tor. A small possible gain asymmetry (|T^ quotes 0.5%) for different zenith 



angle bins can result is a somewhat amplified difference in rates because the 
energy spectrum is rather steep close to the threshold energy (6.5 MeV). The 
energy calibration was done using electron LINAC, which at that time could 
shoot electrons only downwards and hence could not study the asymmetry 
well enough. The gain asymmetry is known to exist from the study of decay 



electrons in the cosmic ray muon data p3 ^s well as in spallation events 



|28[.p| We assume that this will not be an important systematic effect for 
Borexino or KamLAND because the energy deposit is basically isotropic (no 
directional capability) and hence the asymmetry in the detector should not 
result in a systematic effect in the day-night asymmetry. 

The next largest systematic effect is the subtraction of background, ±0.2%. 
If the background events are not completely isotropic, the subtraction de- 
pends on the direction and results in a systematic effect. Again at Borexino 
or KamLAND, the lack of directional correlation eliminates this systematic 
effect. 

If we naively drop these two dominant systematic effects, the size of 
the total systematic uncertainty would be less than 0.1%. Of course, the 
sources of background are very different at Borexino or KamLAND. Possible 
differences in the temperature or Rn level between the day and night times 
could introduce new systematic effects, while our analysis assumed the same 
background level for day and night. This difference, however, can in principle 
be measured using the Bi-Po coincidence. Spallation background (such as 
^^C) should not change between day and night. 

Additionally, the experiments will need to consider other effects, such 
as the contribution of other neutrino sources or the uncertainty in the elec- 
tron number density profile of the Earth. (More on the latter in the next 
section.) We also did not include in our analysis the contribution of neu- 
trinos produced in the CNO cycle, which is about 10% of that from the 
^Be neutrinos. Although we cannot accurately predict the total systematic 
uncertainty at Borexino or KamLAND, we nonetheless find it encouraging 
that the dominant uncertainties at SuperKamiokande are unlikely to affect 
these experiments. 



^The gain asymmetry is now accurately measured using the -'^^N source calibration and 
will be reduced dramatically p3| . 
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4 Measuring the Oscillation Parameters 

In this section, we discuss the possibihty of measuring the value of Am^, sin^ 6 
in the advent of a large day-night effect. In order to do this, data was 
simulated for Am^ = 1.12 x 10^^ eV^, siv? 6 = 0.398, which is close to the 
LOW MSW solution to the solar neutrino puzzle p. For a plot of the "data" 
with different binning options, see Fig. ^. 

In order to deal with the SSM solar neutrino flux and the background 
event rate, we will conservatively "measure" both the background rate and 
the incoming solar neutrino flux by analyzing the seasonal variation [|1^ of 
the day-time data only. This measurement procedure will be incorporated 
in a four parameter x^ analysis (the parameters are Am^, sin^ ^, the solar 
neutrino flux s, and the background rate b) of the data. Explicitly, 

, , jv (data^^^-theof)' ^^ (dataf ^ - theoH ' 



datafsi^M J=i Udataf^ 

(4.1) 
where data"'^ is the night-time "data" binned into N night bins (as de- 
scribed in Sec. |^), data^'^-*' is the day-time "data" binned into M "seasonal 
bins" {e.g. j = 1,2, ...12 months) as described in [|12|. theof'^ is the 
prediction for the number of evens in the i-th night bin. 



theo?'' = 365 



6+.(p°r+(i-p"r)/2) /., (4.2) 



similar to Eq. (|3.1| ). b is the background rate in events per day and s is 
the number of events per day induced solar neutrinos according to the SSM 
prediction for the solar neutrino flux. Similarly, theo^*^*^ is the prediction for 
the day-time flux in the j-th seasonal bin (see W^), 



theof"^ 




^ ^ „ Pee + (1 - Pee)R ' 

(1 — ecos(27rt/year))2 



9j, (4.3) 



where Qj is the number of days in the j-th bin and e = 0.017 is the eccentricity 
of the Earth's orbit. 

It is simple to minimize x^ with respect to s and b, given that x^i^^ b) is a 
quadratic function. The minimization with respect to Am^ and sin^ 9 is done 
numerically. Fig. ^depicts the extracted contours in the (Am^, sin^ 6')-plane, 
in the case of 1 night bin and 10 "uniform" night bins, respectively. 
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Figure 6: Measured values of (A?n^, sin^ 9) at KamLAND after three years of 
running. The data was generated for Am^ = 1.12 x 10^^ eV^, sin^ 9 = 0.398 
(marked with the "star"). The regions obtained by using one night bin and 
ten uniform night bins are shown. 

As Fig. I demonstrates, in the case of 1 night bin, one extracts values of 
Am^ and sin^ 9 which fall into "rings" which correspond roughly to Adn = 
y4^^ ± AA^^ff, where ^4^^^ is the value of the day-night asymmetry for the 
input value of Am^, sin^ 9. In the case of more than one uniform bin, the 
ring degeneracy is broken, and a much more precise determination of the 
oscillation parameters is possible. This is expected, since for Am^ in this 
range the regeneration effect in the Earth exhibits a strong zenith angle 
dependence, as one can easily verified by looking at Fig. ^. 

It is important to note that in the above analysis only statistical uncer- 
tainties were included, while in a real experiment one definitely will have 
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to account for systematic effects as well. In particular, one will need to 
address the uncertainty in the Earth model used in the fit. In producing 
Fig. ^ the same Earth model |2^ was used in generating the "data" and 



in the fit procedure. To understand the effect of using a "wrong" Earth 
model, we have repeated the above analysis using different Earth models in 
the fit. We found the results very encouraging. Even in the case when we 
used for the Earth profile a crude two-step model (a uniform density in the 
mantle and a uniform density in the core), the minimum of x^ occurred at 
lS.rv? = 2.5 X 10~^ eV^, sin^ ^ = 0.24, not far away from the true (input) 
value. Moreover, the x^ value at the minimum was much larger than the 
case with the "true" model (Ax^ = 183 for 18 d.o.f.). This means that in 
a real experiment one will be able to adjust the Earth's model to achieve 
a better fit to the data. Because of the steep rise in x^ value as the Earth 
model is varied, the resulting x^ contours in the (Am^,sin^^) parameter 
space should not be significantly larger than the ones presented here, where 
the Earth model is not varied. As a byproduct of the measurement of the 
neutrino oscillation parameters, it might be possible to use the regeneration 
data to study the interior of the Earth! 

5 Conclusions 

We have studied the effect of the Earth matter on ^Be solar neutrinos. We 
made use of an enlarged parameter space < ^ < 7r/2 and presented 
the sensitivity reach of the KamLAND and Borexino experiments in this 
space. Our results show that both experiments will be sensitive to the Earth 
regeneration effect in a large region which extends into the traditionally 
neglected 9 > 7r/4 part of the parameter space. In particular, for the LOW 
solution one expects to see a greater than 5a effect. On the other hand, both 
experiments will see no day-night effect for the SMA solution and virtually 
no effect for the LMA solution. 

If the experiments see a large Earth regeneration effect, it will be a 
powerful "smoking gun" signature of neutrino oscillations. Furthermore, as 
we have demonstrated, the results of the experiments can be used to measure 
the oscillation parameters. By studying the full zenith angle distribution, 
rather than the usual day-night asymmetry information, one might be able, 
in the case of the LOW solution, to perform a spectacular measurement of 
the parameters. In addition, it might be possible to use the zenith angle 



21 



information to learn about the Earth electron density profile. 

If, on the other hand, no Earth regeneration effect is detected, by com- 
bining this information with the flux measurement from seasonal variation of 
the event rate [0 , a large portion of the parameter space can be excluded. 
If the measured value of the ^Be neutrino flux is large, the exclusion will be 
independent of a specific solar model. 

Both the measurement of the oscillation parameters and the exclusion 
will require a thorough understanding of the systematic uncertainties. We 
have commented on some possible sources of such uncertainties in this paper. 

Overall, Borexino and KamLAND will provide crucial information about 
the solar neutrinos. Not only will the experiments measure the flux of the 
^Be solar neutrinos, but they will also be able to establish or exclude, without 
relying on solar models, the LOW solution based on the Earth regeneration 
effect and the vacuum oscillation solution based on the observed seasonal 
variation of the event rate. Together with results from SuperKamiokande, 
SNO, and the KamLAND reactor neutrino experiment, this information can 
be used to finally unravel the 30-year-old solar neutrino puzzle. 
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A Matter Oscillations and No Level Crossing 

In this appendix, we discuss the survival probability of solar electron neu- 
trinos outside the Sun, in particular the case of no level crossing, i.e., when 
Am^ cos 2^ < in the language of the two neutrino mixing scenario. 
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In the literature, matter effects in the Sun are always considered when 
there is "level crossing" inside of the Sun, i.e., when the light neutrino is pre- 
dominantly of the electron type, and, due to neutrino-electron interactions, 
when the instantaneous Hamiltonian eigenstate with the largest eigenvalue is 
predominantly of the electron type in the Sun's core. The other case, when 
the /ieaf?/ neutrino is predominantly of the electron type, has not been studied 
in the literature in the context of two neutrino oscillations. The authors of 



T6| , however, have considered this possibility in the context of three- flavor 
oscillations. 

The reason for this apparent neglect is simple, and will become clear as 
our results are presented. What happens is that, in the case of no level 
crossing, the electron neutrino survival probability is always bigger than 1/2, 
and therefore it seems that this scenario is not relevant to the solar neutrino 
puzzle. This, however, may not be the case [^. 

Before presenting the expressions for the electron neutrino survival prob- 
abihty outside of the Sun, it is necessary to clearly define the neutrino 
eigenstates and mixing angles. The neutrino mass eigenstates are defined 
in Eq. ( p.l|) , and the notation introduced in Sec. ^ will be used. In what 
follows our convention is Arn^ > and < sin^ 9 < 1. 

Inside the Sun the Hamiltonian has the form 

- V2GFNn{r)/2 

^cos2e + V2G,N^{r) ^sin20\ 

^sm29 4^ cos 2^ ) ' ^ ^^ 

were Pi, is the solar neutrino momentum, Ne{r) {Nn{r)) is the electron (neu- 
tron) number density at a distance r from the Sun's core, m"^^^ = ml + ml 
and Am^ = m,\ — mj\, m| and m^ being the mass eigenvalues in vacuum. The 
eigenvalues of the instantaneous Hamiltonian are 




^^' ~ ^"^ Ap, ^ 2 2 

+ [V2GFN^ir)y - ^^V2GFNeir) cos 2^ 



fAm^ 

1/2 



(A.2) 



For the study of neutrino oscillations, terms common to both states are irrel- 
evant, and the first three terms can be dropped. One can also safely replace 
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Pu by El, in the remainder. The instantaneous Hamiltonian eigenstates in 
terms of flavor eigenstates are 

|z/_(r)) = cos 9 M{r)\ue) - sm9M{r)\Uf,), (A. 3) 

|z/+(r)) = sin^M(r)|z/e) + cos^M(r)|z/^). (A.4) 

Here OMi^) is the matter mixing angle, given by 

cos2^m(^) = 

^m'^cos2e-2E^y/2GFNe{r) 

^(Am2)2 + {2E^V2GFNe{r)Y - AAm'^E^V^GpNeir) cos 2^' 

Assuming that A^e -^ as r — i> oo, it is easy to see that \i'+{r -^ oo)) -^ |z/2) 
and |z/_(r -^ oo)) -^ |z/i). Therefore, if the transition from the core of the 
Sun to vacuum is adiabatic, a state which is created as a |z/+(0)) (|z^-(0))) 
will exit the Sun as a |z/2) (l^i))- 

Having established the notation, it is very easy to estimate the survival 
probability for electron neutrinos that are created in the Sun's core and are 
detected on the Earth, in the limit that Am'^/2Eiy is much smaller than the 
Earth-Sun distance, such that oscillations in vacuum between ui and z/2 states 
are "averaged out." There are four possible "propagation paths" that the 
solar neutrino can follow: 

z/e -^ z/+(p = sin^ 9m) -^ i^2{p = '^ - Pc) ^ t^e{p = siu^ 0) 

or 
z/e -* u^{p = sin^ 6m) —^vi{p = Pc) -^ Veip = cos^ 6) 

or (A.6) 

Z/g -* J^'iP = COS^ 6m) ^ l^lip = 1 - Pc) ^ l^eip = COS^ 6) 

or 

Ue -^ '^'{P = COS^ 6m) -^ 1'2{P = Pc) -^ l^eip = siu^ 6), 

where p is the probability that a given "step" takes place, 6m = 6'm(0) and Pc 
is the jumping probability, i.e. the probability that during the evolution from 
the Sun's core to vacuum the neutrino changes from one set of instantaneous 
Hamiltonian eigenstates to the other. 

Therefore, the probabilities of finding the mass eigenstates z/i and z/2 far 
from the Sun are given by 

Pi = sm^6MPc + cos^6M{l-Pc) (A.7) 

P2 = I -Pi, (A.8) 
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where 9m is that at the production point, [J] and the electron neutrino survival 
probability (-Pee) at the surface of the Earth is 



Pe, 



= Pi cos^ + P2 sin^ ^ 

= sin^ ^m((1 - Pc) sin^ 6 + P^ cos^ 6) 

+ cos^ ^m((1 - Pc) cos^ 9 + Pc sin^ 9). 



(A.9) 



All equalities hold as long as two mass eigenstates appear as an incoherent 
mixture (true for Am^ ^ 10^^ eV^ for ^Be neutrinos). In deriving Eq. ( A.9 ), 
no assumption was made with respect to the value of cos 2^, and therefore it 
should be valid for the entire range of < sin^ 9 < 1. 

Given Eq. ( |A.9[ ), it is easy to show that for 9 < ii/A, P^e can be (much) 
smaller than 1/2, while for 9 > tt/4, P^e is larger than 1/2 (indeed, it will be 
shown that Pee > P^^,, the (averaged) vacuum survival probability). 

First, note that —1 < cos29m < cos 2^. The equalities are saturated 
when V^GpNeiO) > Am'^/2E^ or V^GpNeiO) < Am^/2E^, respectively. 
More quantitatively 



Am^ 



2E,V2GfNM 



0.98 



Am^ \ 
lO-^eV^j 



0.862MeV\ 



E,. 



-) 



for an average core electron number density of 79 moles/cm'^ 
in the case of ^Be neutrinos and Am'^ -C 10~^ eV^, 



(A.IO) 
rU| . Therefore, 



cos 2^ 



M 



1 / Am2sin26' 

2 [2E,V2GFNe{0), 



O 



Am? 



,2E,V2GfN,{0), 



and 



P ~ 

-* ee — 



(I-Pjsin^^ + Pecos^e. 



We will soon show that Pc G [0, cos^ 6'],Q so that, in the limit cos 29 

Pee e [sin^ 9, sin^ 9 + cos^ 9] or 
Pee G [sin^ 9 + cos^ 9, sin^ 9] 



M 



(A.ll) 
(A.12) 

(A.13) 
(A.14) 



"in our numerical analyses, we integrate over the production region using the profile 
given in [p2| . The interference between 1/+ and j/_ states in Eq. ( A. 6 ) vanishes upon 
averaging over the production region independent of Arri^ or energy. 

**This is not hard to see. It is known that, if Arri^ is large enough, the adiabatic 
approximation should hold, and therefore Pc ^ for large enough Arn^. On the other 
hand, if Arri^ is small enough, one should reproduce the vacuum oscillation result (as 



in the just-so scenario), and, from Eq. (A.9), it is easy to see that this happens when 
Pc — > cos-^ 6 and cos26'm ^ — 1. 
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Eq. ( |A:T^ ) (Eq. i^^J§) applies if sin^^ < cos^ ^ (sin^ ^ > cos^^). This is 
easy to see because sin'^ 9 + cos^ 9 = 1 — (1/2) sin^ 29 is the vacuum survival 
probability P^^ and 

P;, = l-2sin2^(l-sin2 0) 

= l-2sin2e + 2sin^^, (A. 15) 

which is bigger (smaller) than sin^ 9 if sin^ 9 < cos^ 9 (sin^ 9 > cos^ 9). 

When \/2GfN(,(0) <^ /S.m? /2Ejj matter interactions should be irrelevant, 
and it is easy to see from Eq. ( |A.5| ) that cos29m -^ cos 26'. In this limit 
Pc — ^ 0, since we are deep into the adiabatic region (as will be shown later) 
and Pee —>■ Pee- Bcforc summarizing the behavior of the ^Be electron neutrino 
survival probability we will determine some expression for Pc. 

Assuming an exponential profile for the electron number density inside 
the Sun (A^e('") = A'e(O) exp(— r/ro)), Schrodinger's equation can be solved 
analytically pT| , pT[] , and it is shown that, in the range of the neutrino 



oscillation parameter space relevant for addressing the solar neutrino puzzle, 
Eq. (|A.9|) is indeed a very good approximation for Pee and that Pc is given 
by|2T],|| 

P^= i_,-7 - (A-16) 

where 

^ Am2 ^^^/ Am2 \ /0.862 MeV\ ,,^^, 

for ro = Po/10.54 = 6.60 x 10^ km. 

According to the author of [^, Eq. (|A.16|) only holds for Am^ cos 29 > 0. 



We will prove shortly, however, that Eq. ( [A.16| ) also applies in the case of 
no level crossing, when the heavy mass eigenstate is predominantly of the 
electron type, i.e. when sin^ 6' > cos^ 6'. Assuming that this is indeed the 
case, we can finish our discussion on the behavior of the electron neutrino 
survival probability, using ^Be neutrinos as an example. 

When Am^ -C 10~^ eV^, cos26'm = —1 and Pc = cos"^ 9. In this case 
we argued and one can explicitly check that Pee = Pe^'e-EH For 10~^ eV^ -C 

't^^Indeed, this is the region of the "just-so" solution. As a matter of fact, in this region 
the distance dependent vacuum oscillations do not average out when the neutrinos are 
detected at the Earth, and one should use the appropriate position dependent expression. 
It is known, however, that the equality Pee{r) = Pee^'"') holds, up to a phase |l2[ ^. That 
this is also true for cos 26* < was explicitly checked starting with the exact solutions to 
Schrodinger's equation [|T[ pi[ . 
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Am^ < 10~^ eV^, cos2^m = -1 and Pc -^ 0. In this case Pee — sin^^. 
This is the adiabatic region. For Arn^ ^ 10^^ eV^, matter effects become 
irrelevant and cos 20^/ = cos 26*, Pc = 0. Again Pee = Pee- Therefore, 
Eqs. ( |A.13| , [ATT4D apply for all values of interest, and one can get a very large 
suppression of Pee if sin^6' -C 1. On the other hand, in the case of no level 
crossing. Pee is always bigger than PJ^ > 1/2. 

Fig. depicts the behavior of Pg^e^*^-* as a function of Am^, for different 
values of the vacuum mixing angle. The preferred values from the overall 
rate analysis at the Homestake, Sage and Gallex, and SuperKamiokande 
experiments are indicated by stars. The four plots are labeled SMA, 
LMA, LOW to indicate that they contain the best fit values of 6 for the Small 
Mixing Angle, Large Mixing Angle and LOW Am"^ solutions [^], respectively, 
and INT to indicate an intermediate value of 6 between the SMA and LMA 
solutions. The dotted line indicates the value of P^e- Similarly, Fig. |^ depicts 
Pje^'^^ as a function of sin^ 9 for different values of the mass squared difference. 
We use the same notation as the one used in Fig. |^, and the vertical dashed 
lines indicate the mixing angle for maximal vacuum mixing (sin^ 6' = 1/2). 
Note that at this point P£^^) = PJ^ = 1/2. 

Finally, we argue that Eq. ( |A.16| ) holds for all values of cos 2^. When 
\/2GFNe{0) ^ |Am^|/2_E'j,, it is very simple to derive Pee, following the 
exact solution |^, ^ to Schrodinger's equation and taking the appropriate 



limits. According to Eq. (39) in [^ 
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Figure 7: The electron neutrino survival probability as a function Am^, for 
different values of the vacuum mixing angle, namely, cos 2^ = ±0.997 (SMA), 
cos 2^ = ±0.8 (INT), cos 20 = ±0.58 (LMA), and cos 20 = ±0.24 (LOW). 
The upper (lower) lines are for the negative (positive) sign of cos 26* < 0. 
The stars indicate the preferred points from the overall rate analysis of the 
existing data , and the horizontal dotted lines indicate the vacuum survival 
probability, Pl^ = 1/2 - 1/2 sin^ 20 . 
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Figure 8: The electron neutrino survival probability as a function sin^ ^, for 
different values of Am^, namely, Am^ = 1.3 x 10"^ eV^ (LOW), Am^ = 5.0 x 
10-6 eV2 (SMA), Am^ = 1.4x10-^ eV2(LMA), and Am^ = 1x10"^ eV^. The 
stars indicate the preferred points from the overall rate analysis of the existing 
data 0, and the dashed lines indicate the vacuum survival probabihty, P^^ = 
1/2 - 1/2 sin^ 29. 
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where ho = Am? /2Ey and Pc is given exactly by Eq. (A.l(j). Therefore 



Pee = l-P^e^{l- Pc) sin^ e + P, COS^ 9. (A. 19) 



Since, in deriving Eq. ( [A. 181 ), no assumptions with respect to the sign of Am 



2 



or the value of 6 were made, it should be applicable in all casesj^ as long as 
V2GFNe{0) > \Am'^\/2E^. Indeed, from Eqs. O ^) it is easy to note 
that, in the limit V^GpNeiO) > \Am'^\/2E^, cos29m -^ -1 and Eq. dTTTOD 
is exactly reproduced. 

It should be noted that, in the region of no level crossing, Pc ~ for 
the entire Am'^ range that we are interested in. Indeed, only in the region 
of the "just-so" solution the crossing probability is appreciably different 
from zero, and it was in this region (cos26'a/ = —1) that we explicitly 
showed (Eq. ( |A.18| )) that Eq. ( |A.16| ) holds. If there are any interesting 
phenomenological consequences for the 6 > tt/4 case in the "just-so" to 
MSW transition region of the parameter space still remains to be seen [^. 
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